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ABSTRACT 


The  goal  of  this  project  was  to  improve  our  understanding  of  electromagnetic  scattering  from 
conducting  rough  surfaces  when  the  roughness  is  small  in  comparison  with  the  radiation  wave¬ 
length.  Research  was  conducted  in  two  related  areas  for  p-polarized  scattering  from  conducting 
surfaces:  (1)  the  generation  of  surface  waves,  and  (2)  the  development  of  renormalized  perturba¬ 
tion  theory  for  scattering  from  perfectly  conducting  surfaces.  To  understand  the  excitation  of  sur¬ 
face  waves  better,  a  numerical  method  was  developed  for  displaying  the  field  near  the  surface; 
surface  waves  can  then  be  shown  explicitly.  The  method  is  based  on  integral  equation  solutions, 
and  no  fundamental  approximations  are  required.  The  existence  of  surface  waves  for  p-polarized 
scattering  leads  to  difficulties  in  theoretical  treatments  of  scattering  from  conductors.  This  work 
is  specialized  to  the  case  of  scattering  from  perfect  conductors,  for  which  the  theoretical  problems 
are  most  evident.  A  renormalized  perturbation  theory  is  implemented  that  is  free  from  the  defi¬ 
ciencies  inherent  in  standard  perturbation  theory.  The  renormalized  perturbation  theory  is 
extended  beyond  lowest  order,  and  comparisons  with  exact  integral  equation  results  show  that  this 
approach  yields  predictions  for  the  bistatic  scattering  cross  section  that  are  accurate  over  a  broad 
range  of  surface  parameters. 
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1.  Problem  Statement 

The  goal  of  this  project  was  to  improve  our  understanding  of  electromagnetic  scattering  from 
conducting  rough  surfaces,  when  the  roughness  is  small  in  comparison  with  the  radiation  wave¬ 
length.  Research  has  been  conducted  on  two  related  areas  for  p-polarized  scattering  from  con¬ 
ducting  surfaces:  (1)  the  generation  of  surface  waves,  and  (2)  the  development  of  renormalized 
perturbation  theory  for  scattering  from  perfectly  conducting  surfaces.  For  both  topics,  use  has 
been  made  of  the  exact  integral  equation  method  restricted  to  one-dimensional  (1-D)  surfaces,  and 
the  restriction  to  1-D  surfaces  has  been  made  throughout  in  the  analytic  development  as  well. 

Surface  waves  can  be  generated  when  p-polarized  light  (electric  vector  in  the  plane  of  inci¬ 
dence)  scatters  from  good  or  perfect  conductors  with  small  roughness.  To  understand  the  excita¬ 
tion  of  surface  waves  better,  a  numerical  method  has  been  developed  for  displaying  the  field  near 
the  surface;  surface  waves  can  then  be  shown  explicitly.  The  method  is  based  on  integral  equation 
solutions,  and  no  fundamental  approximations  are  required.  An  outline  of  the  numerical  method 
and  a  discussion  of  results  will  be  given. 

The  existence  of  surface  waves  for  p-polarized  scattering  leads  to  difficulties  in  theoretical 
treatments  of  scattering  from  conductors.  Here  we  specialize  to  the  case  of  scattering  from  perfect 
conductors,  for  which  the  theoretical  problems  are  most  evident.  When  standard  perturbation  the¬ 
ory  is  applied  to  this  case  (the  Neumann  boundary  condition  problem),  lowest-order  results  have 
unphysical  properties  near  grazing  (i.e.,  at  angles  of  incidence  or  scattering  near  ±90°)  and  there¬ 
fore  must  be  incorrect  near  grazing.  When  standard  perturbation  theory  is  taken  to  second  order, 
the  expression  for  the  scattering  cross  section  diverges.1  Building  on  theoretical  work  of  A. A. 
Maradudin  and  colleagues,2'5  a  renormalized  perturbation  theory  has  been  developed  and  imple¬ 
mented  that  is  free  from  the  deficiencies  just  mentioned.  The  goal  has  been  to  develop  a  perturba¬ 
tion  theory  with  an  implementation  capable  of  giving  accurate  results  for  kh  up  to  0.5,  where  k  is 
the  radiation  wavenumber  and  where  h  is  the  rms  height  of  the  surface  roughness.  Perturbation 
theory,  if  taken  beyond  lowest  order,  is  valid  for  s-polarized  scattering  from  perfect  conductors 
(Dirichlet  boundary  condition)  up  to  about  this  level  of  roughness.6  Comparisons  with  exact  inte¬ 
gral  equation  results  show  that  this  goal  has  now  been  essentially  reached  for  the  Neumann 
boundary  condition.  An  outline  will  be  given  of  the  renormalized  perturbation  theory  followed  by 
a  discussion  of  implementation  issues  and  a  comparison  with  integral  equation  results. 


2.  Field  Visualization  in  P-Polarized  Scattering 

We  consider  scattering  from  both  perfect  and  good  conductors.  For  perfect  conductors,  sur¬ 
face  roughness  is  required  for  surface  waves  to  exist.  For  good  conductors,  surface  roughness  can 
couple  energy  from  volume  propagating  waves  into  surface  waves,  but  after  the  surface  wave  has 
been  generated,  roughness  is  not  necessary  for  the  surface  wave  to  exist.  By  displaying  the 
near-surface  fields,  we  can  study  the  properties  of  surface  waves  for  both  good  and  perfect  con¬ 
ductors. 

2.1  Numerical  method 

n 

The  integral  equations  given  below  are  equivalent  to  those  derived  by  Maradudin  et  al.  We 
are  concerned  with  the  scattering  of  p-polarized  light  from  a  rough  metal  surface.  Time  harmonic 
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fields  are  assumed,  and  the  time  dependence  is  omitted.  The  plane  of  incidence  is  the  x-z  plane, 
and  the  rough  surface  profile  is  defined  by  z  =  fix),  where  the  mean  surface  height  is  zero  and 
where  the  z  axis  is  directed  upward  from  the  metal  into  a  vacuum.  For  p-polarized  scattering  from 
a  1-D  surface,  the  y  component  is  the  only  nonvanishing  component  of  the  magnetic  field,  which 
is  denoted  by  H( r)  where  r  is  a  2-D  position  vector.  We  wish  to  solve  for  H(r)  in  the  vicinity  of 
the  surface  when  a  field  Hmc(r)  is  incident  on  the  metal  from  the  vacuum  side  of  the  interface. 


If  the  metal  is  a  perfect  conductor,  the  field  vanishes  inside  the  conductor,  and  the  boundary 
condition  on  the  rough  surface  is 


r)  _ 

dn 


(1) 


where  d/dn  =  n  V,  with  n  a  unit  vector  normal  to  the  surface  and  directed  from  the  metal  to  the 
vacuum.  The  solution  of  the  following  integral  equation  gives  H( r)  on  z  =  fix) 

r3Gn(r,  r')  #  , 

H( r)  =  2H.nc(r)  +  2J— ^7 - H( r')  ds'  •  (2) 

S 

In  (2)  both  r  and  r'  are  confined  to  the  rough  surface,  and  ds'  =  y(x')dx',  where  ^(x)  =  1  +  (df(x)/ 
dxj2.  Also  in  (2)  G0  is  the  free  space  Green’s  function  in  the  vacuum 

G0(r,r')  =  ij^Vlf-r'l)  ,  (3) 


where  H0(^  is  the  zero-order  Hankel  function  of  the  first  kind  and  where  k  =  to/c.  The  incident  field 
is  taken  as  a  Gaussian  tapered  plane  wave.8  The  scattered  field  at  arbitrary  points  above  the  surface 
is  then  given  by 


*,(r) 


,8Gn(r,  r') 

=  |  °V>— 


(4) 


In  scattering  cross  section  calculations,  the  scattered  field  is  needed  in  the  far  field,  and  a  far  field 
form  of  (4)  is  used.  In  order  to  obtain  the  scattered  field  close  to  the  surface,  however,  the  more 
general  form  given  by  (4)  is  used  here. 


If  the  metal  is  a  good  conductor,  the  magnetic  field  will  penetrate  into  the  conductor,  and  we 
must  deal  with  the  more  general  two-medium  problem  separated  by  a  rough  interface.  The  bound¬ 
ary  conditions  at  the  interface  are  that  H{ r)  and  ( 1/e)  [dH  (r)  /dn]  are  continuous;  £  is  the 
dielectric  constant,  which  here  is  unity  above  the  surface.  The  following  coupled  integral  equa¬ 
tions  are  obtained  for  H{ r)  and  its  normal  derivative 


»(r>  =  2Hjnc(r)+2| 


dG0,r’  r’>  r  ,,  3W(r') 

dn’  ff(r)~G0<r-r)-Srj 


ds' 


(5) 


and 
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H(r)  =  -2j 


3G2(r,r')  ,,3H(r') 

1  H{r')  -e2G2(r,r  ) — ^ 


dn' 


dn' 


(6) 


In  (6)  e2  is  the  complex  dielectric  constant  of  the  metal  and 

G2(r,r')  =jH0(1)<*2|r-r1)  (7) 

where  k2  =  Je~2  oo/c  .  Note  that  (6)  requires  the  evaluation  of  Hankel  functions  of  complex 
argument.  After  solving  (5)  and  (6)  for  the  surface  fields,  the  scattered  field  in  the  vacuum  above 
the  surface  is  given  by 


Hs(r) 


dG0(r,  r') 
dn' 


H(r')-GJr,  r') 


<m  O 

dn' 


ds' 


’  (8) 


The  integral  equation,  (2),  or  the  coupled  integral  equation  system,  (5)  and  (6),  can  be 
reduced  to  matrix  equations  and  solved  numerically  following  standard  methods  (see,  for  example 
Refs.  7  and  8).  The  solution  of  (2)  gives  H( r)  on  the  surface  z  =/(*)  at  a  set  of  N  equally  spaced 
points  along  the  x-axis.  For  the  examples  to  be  discussed,  we  have  used  an  x-axis  partition  interval 
Ax  =  A/ 10,  where  X  is  the  wavelength  in  vacuum.  Equations  (4)  and  (8)  then  reduce  to  sums  over 
these  N  points.  Either  (4)  or  (8)  is  used  to  obtain  the  scattered  field  at  a  rectangular  array  of  field 
points  close  to  the  surface  for  display  in  color-coded  plots. 

Two  important  modifications  were  made  to  the  general  procedure  just  described.  The  first  is 
necessary  only  for  very  good  conductors  of  the  type  to  be  considered  here,  but  the  second  is  nec¬ 
essary  for  all  cases. 

Scattering  from  gold  surfaces  at  a  wavelength  of  3.392  pm  has  been  chosen  to  illustrate  sur¬ 
face  wave  properties.  Gold  is  a  very  good  conductor  for  radiation  of  this  wavelength,  and  yet,  as 
we  will  see,  the  surface  wave  properties  for  gold  are  very  different  than  those  for  a  perfect  con¬ 
ductor.  The  complex  index  of  refraction  for  gold  at  3.392  pm  is  n  =  1.958  +  20.7 /,  and  the  corre¬ 
sponding  complex  dielectric  constant  is  £2  =  —424.7  +  81.1/.  For  such  values  of  e2  it  is  not 
difficult  to  show  that  G2  will  vary  significantly  over  Ax  even  when  Ax  =  A/10.  Thus,  the  normal 
practice  of  approximating  G2  and  its  normal  derivative  as  constant  over  each  Ax  when  developing 
the  matrix  equations  leads  to  significant  error.  Reducing  Ax  to  X/40  gives  improved  results;  how¬ 
ever,  even  finer  partitioning  appears  necessary  for  this  case,  and  the  approach  of  refining  the  par¬ 
titioning  becomes  computationally  unwieldy.  One  approach  that  does  work  in  evaluating  the 
matrix  elements  is  to  approximate  H( r)  and  dH(  r)  /dn  in  (6)  as  constant  over  Ax,  and  then  to 
integrate  the  more  rapidly  varying  quantities  G2  and  dG2/dn  numerically  over  Ax.  This  method 
has  been  developed  and  gives  reliable  results,  but  a  simpler  method  to  be  described  next  gives 
essentially  identical  results  and  therefore  has  been  used  for  the  examples  to  be  discussed. 
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The  conductivity  of  gold  at  3.392  pm  is  sufficiently  high  that  H( r)  and  its  normal  derivative 
are  very  accurately  related  by  an  impedance  condition 

(9) 

The  impedance  condition  (9)  is  just  the  first  term  in  an  expansion  in  powers  of  skin  depth,9 
but  it  is  entirely  adequate  for  the  present  case.  Use  of  (9)  in  (5)  leads  to  a  single  integral  equation 
for  H( r),  and  (6)  is  not  used.  Thus,  the  matrix  size  for  our  system  of  equations  is  reduced  from  2N 
x  2 N  to  Nx  N,  which  yields  a  major  reduction  in  computation  time.  In  addition,  only  Hankel 
functions  of  real  arguments  are  required,  which  is  an  added  benefit.  After  H( r)  is  obtained  from 
(5),  the  scattered  field  is  found  with  (8),  again  using  (9).  Comparisons  between  the  solutions 
obtained  using  the  full  coupled  equation  system  of  (5)  and  (6)  and  using  the  impedance  condition 
(9)  in  (5)  show  essentially  perfect  agreement  for  gold  at  3.392  pm.  Therefore,  the  impedance 
method  was  used  for  the  good  conductor  cases  to  be  discussed.  For  other  metals  or  frequencies, 
such  agreement  would  need  to  be  reconfirmed,  and  an  improved  impedance  approximation  or  the 
full  solution  method  may  well  be  necessary. 

The  second  modification  to  the  general  procedure  for  calculating  the  scattered  field  is  made 
only  for  field  points  very  close  to  the  surface.  Typically,  (4)  or  (8)  is  used  to  find  Hs(r)  in  a  rectan¬ 
gular  array  of  field  points  with  horizontal  and  vertical  spacings  of  A/ 10.  The  procedure  as 
described  works  well  except  for  field  points  that  happen  to  be  within  a  small  fraction  of  a  wave¬ 
length  from  one  of  the  set  of  N  points  that  define  the  rough  surface.  The  field  is  compute  on  a  pre¬ 
scribed  set  of  points  independent  of  the  rough  surface  profile,  and  when  a  field  point  is  very  close 
to  a  surface  point,  the  approximation  of  replacing  (4)  or  (8)  as  a  simple  sum  becomes  inaccurate. 
The  rapid  variation  of  G0  as  r'  passes  close  to  r  must  be  taken  into  account  in  the  integration  to 
obtain  accurate  results. 

The  modified  procedure  is  as  follows.  The  medium  in  which  each  field  point  resides  is  deter¬ 
mined  assuming  that  straight  lines  connect  the  N  points  defining  the  surface.  If  the  field  point  lies 
inside  the  metal,  the  field  is  not  computed,  and  the  color  set  to  black.  For  the  good  conductor  case, 
the  actual  field  can  be  found  using  an  expression  analogous  to  (8)  that  is  valid  in  the  second 
medium.  However,  for  gold  the  field  drops  off  very  rapidly  with  increased  depth  in  the  metal,  and 
for  simplicity  this  field  has  not  been  included.  Next,  the  field  point  is  tested  for  closeness  to  sur¬ 
face  points.  Whenever  k\  r  -  r'l  <  0.9,  the  corresponding  term  in  (4)  or  (8)  is  replaced  by  an  ana¬ 
lytic  integration  over  Ax.  This  is  done  by  making  a  small  argument  expansion  of  HQ^  and 
keeping  terms  up  to  the  argument-squared  term  (the  first  four  terms).  The  integration  over  AlX  can 
then  be  done  analytically,  taking  the  normal  derivative  of  H{ r)  as  constant  over  this  interval.  Sim¬ 
ilarly,  the  normal  derivative  of  G0  yields  Hyl)  which  is  also  expanded  in  a  small  argument  series 
with  terms  kept  up  to  the  argument  squared  term  (the  first  five  terms).  Again,  an  analytic  integra¬ 
tion  is  done  over  Ax.  The  result  of  this  modification  is  to  produce  a  smooth  scattered  field  right  up 
to  the  surface,  including  the  regions  within  the  corrugations  of  the  rough  surface. 
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2.2  Examples  of  field  visualization  in  p-polarized  scattering 

Examples  of  fields  near  rough  surfaces  for  good  and  perfect  conductors  are  shown  in 
Figs.  1-3.  These  color  field  maps  were  made  by  first  contouring  the  field  values  that  were  com¬ 
puted  as  described  in  Section  2.1  and  then  assigning  colors  to  the  field  values  on  a  linear  scale. 
For  the  examples  shown,  the  surface  was  modeled  with  a  Gaussian  surface  height  spectrum  with 
kh  =  0. 1  or  0.5  and  kl  =  1.0,  where  k  is  the  radiation  wavenumber,  h  is  the  rms  surface  height,  and 
/  is  the  He  correlation  length.  The  rough  surface  extends  from  -450  to  +450  on  the  horizontal 
axis,  and  a  wavelength  equals  ten  units.  The  incident  field  is  a  Gaussian  tapered  plane  wave  cen¬ 
tered  at  zero  on  the  horizontal  scale  and  with  an  incident  angle  of  45°.  The  He  points  for  the  inci¬ 
dent  field  are  at  ±100  on  the  horizontal  scale,  and  here  we  display  the  scattered  field  well  away 
from  the  region  where  the  incident  field  has  significant  magnitude. 

For  kh  =  0.1  most  of  the  energy  in  the  scattered  field  resides  in  the  coherent  component, 
which  propagates  in  the  specular  direction.  The  smaller  incoherent  field  scatters  in  all  directions. 
It  turns  out  that  for  kl  =  1.0,  and  with  the  boundary  conditions  under  study,  more  incoherent 
energy  is  scattered  in  the  back  direction  than  the  forward  direction.  (For  kl  =  2.0,  there  is  about 
equal  backward  and  forward  scattering,  and  as  kl  is  increased  further  forward  scattering  domi¬ 
nates.)  Therefore  we  display  fields  near  the  surface  in  the  back  direction  to  look  for  surface 
waves;  the  fields  in  the  forward  direction  are  similar  but  with  somewhat  reduced  magnitude. 

Figure  1  shows  the  magnitude  of  the  scattered  field  for  a  perfect  conductor  (a)  and  a  good 
conductor  (b)  for  the  same  surface  realization  (kh  =  0.1).  The  roughness  of  the  surface  is  too  small 
to  be  seen  for  this  case,  and  as  mentioned  in  Section  2.1  only  the  field  on  the  vacuum  side  of  the 
interface  is  shown.  In  the  initial  integral  equation  calculation  of  the  surface  field,  the  rough  sur¬ 
face  continues  on  to  the  left  to  -450  to  avoid  edge  effects  in  the  field  maps  shown.  Numerical 
studies  show  that  extending  the  surface  endpoint  farther  to  the  left  makes  no  significant  differ¬ 
ence.  The  incident  field  has  not  been  added  to  the  scattered  field  to  yield  the  total  field  and  thus  is 
not  included.  In  any  case,  the  incident  field  would  be  almost  entirely  out  of  the  field  of  view;  it 
propagates  from  the  upper  left  to  the  lower  right  and  is  concentrated  either  above  or  to  the  right  of 
the  view  shown.  Thus,  the  figures  display  the  near-surface  region  to  the  left  of  where  the  main 
part  of  the  incident  field  illuminates  the  surface.  The  range  on  the  magnitude  scale  is  from  0.0  to 
0.25,  while  the  magnitude  of  the  incident  field  is  1 .0  at  0.0  on  the  horizontal  scale. 

In  Fig.  1  the  large  feature  in  the  middle  to  upper  part  of  both  (a)  and  (b)  is  a  lobe  of  freely 
propagating  scattered  energy  moving  up  and  to  the  left.  This  structure  depends  on  the  particular 
surface  realization  used  and  would  differ  for  each  realization.  Note  that  this  lobe  is  similar  for  the 
good  and  perfect  conductor  models.  In  fact,  evaluation  of  the  scattered  field  in  the  far  field  of  the 
surface  yields  results  that  are  very  close  for  both  models.  However,  the  structure  near  the  surface 
is  very  different  for  these  two  cases,  and  this  difference  can  have  important  consequences  for  the¬ 
oretical  treatments  of  scattering. 

A  surface  wave  should  appear  as  an  enhancement  of  the  field  near  the  surface.  In  Fig.  1  the 
good  conductor  shows  a  clear  surface  wave,  while  the  presence  of  a  surface  wave  is  not  evident 
for  the  perfect  conductor.  Note  also  that  the  field  magnitude  in  the  surface  wave  shows  a  clear 
modulation  with  a  half-wavelength  period  (five  units).  Numerical  studies  show  this  to  be  a  stand¬ 
ing  wave  pattern  set  up  by  right-  and  left-propagating  surface  waves.  Though  a  predominantly 
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Figure  1.  Backscattered  field  magnitude  near  rough  surface  with  kh  =  0.1  and  kl  =  1.0:  (a)  perfect 
conductor,  (b)  good  conductor.  The  view  is  to  the  left  of  the  region  directly  illuminated 
with  a  tapered  plane  wave  incident  field  (0i  =  45°),  which  is  centered  at  zero  on  the 
horizontal  scale.  A  radiation  wavelength  equals  ten  units,  and  the  color  scale  is  linear  in 
the  field  magnitude. 
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left-propagating  surface  wave  might  be  expected  for  the  geometry  of  Fig.  1,  the  roughness  is  suf¬ 
ficient  to  generate  a  right-going  wave  as  well.  If  the  surface  is  set  flat  to  the  left  of -300,  for  exam¬ 
ple,  the  modulation  vanishes  to  the  left  of  that  point  and  grows  rapidly  to  the  right  of  that  point. 
The  modulation  period  remains  at  one  half  wavelength  when  kl  is  increased  to  2.0,  showing  that 
the  scale  length  of  the  roughness  does  not  determine  the  modulation  period.  Analytic  estimates 
show  that  when  an  evanescent  wave  interacts  once  with  the  rough  surface  with  kl  =  1.0,  there  is  a 
greater  amplitude  (or,  more  precisely,  a  greater  transition  matrix  element)  for  scattering  into 
back-traveling  evanescent  waves  than  into  forward-traveling  evanescent  waves.  This  occurs  for 
essentially  the  same  reason  that  the  far  field  backscattering  level  is  higher  than  the  forward  scat¬ 
tering  level.  The  dominance  of  backscattering  apparently  leads  to  the  establishment  of  the  stand¬ 
ing  wave  pattern.  When  kl  is  increased  to  3  or  4,  the  amplitude  for  scattering  into 
forward-traveling  waves,  freely  propagating  or  evanescent,  is  much  greater  than  the  amplitude  for 
scattering  into  backward  traveling  waves.  In  this  case,  surface  waves  are  still  produced  for  the 
good  conductor  (now  looking  in  the  forward  region,  e.g.,  from  +200  to  +400),  but  the 
half-wavelength  modulation  is  absent. 

As  mentioned,  the  perfect  conductor  does  not  show  a  clear  surface  wave  in  Fig.  1(a),  but  on 
close  inspection  there  is  some  evidence  for  one.  There  is  half-wavelength  modulation  of  the  field 
near  the  surface  (and  extending  well  above  it),  and  there  is  some  weak  enhancement  of  the  field 
near  the  surface.  By  studying  a  series  of  cases  with  increasing  kh,  one  can  see  more  clearly  that  a 
surface  wave  does  exist  in  the  perfect  conductor  case,  though  it  is  weak  when  kh  =  0.1  and  not 
particularly  convincing  from  Fig.  1(a)  alone. 

Figure  2  compares  field  magnitudes  for  perfect  [Fig.  2(a)]  and  good  [Fig.  2(b)]  conductors 
when  kh  is  increased  to  0.5  with  kl  and  the  incident  angle  unchanged.  Note  that  the  view  here  is  in 
closer  to  the  surface,  showing  half  as  much  surface  length.  The  roughness  of  the  surface  can  now 
be  seen,  but  the  “blocky”  appearance  of  the  interface  is  an  artifact  of  the  plotting  process  and  does 
not  accurately  represent  the  rough  interface  profile.  (The  potential  problem  here  is  that  the  con¬ 
touring  routine  can  produce  a  blurred  interface.  To  avoid  this,  the  “field  magnitude”  has  been  set 
to  -100  when  the  field  point  lies  in  the  metal.  The  contouring  routine  then  creates  a  sharp  bound¬ 
ary  along  edges  of  the  1  unit  x  1  unit  cells  used  to  sample  the  field.  A  more  accurate  representa¬ 
tion  of  the  interface  can  be  obtained  by  reducing  this  cell  size,  but  a  smaller  grid  is  not  necessary 
to  represent  the  field  itself  accurately.) 

For  the  case  shown  in  Fig.  2,  the  roughness  is  higher  than  for  Fig.  1,  and  consequently  more 
energy  is  scattered  in  the  back  direction.  Because  some  energy  is  lost  in  the  metal  for  the  good 
conductor  (resistive  loss),  the  levels  are  lower.  Note  that  the  field  propagating  at  low  scattered 
grazing  angles,  that  is,  nearly  parallel  to  the  surface,  is  much  reduced  in  magnitude.  Apparently, 
the  field  components  propagating  at  low  grazing  angles  close  to  the  surface  have  a  high  probabil¬ 
ity  of  rescattering  from  the  rough  surface  to  higher  grazing  angles.  The  field  propagating  close  to 
the  surface  is  therefore  effectively  stripped  away. 

It  is  also  evident  in  Fig.  2  that  there  is  field  enhancement  close  to  the  surface  (within  about 
five  units,  or  one-half  wavelength),  and  hence  a  surface  wave  is  present,  though  it  decays  more 
rapidly  with  distance  along  the  surface.  For  this  case,  the  surface  wave  field  structure  is  similar  for 
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the  good  and  perfect  conductor  models.  Again,  a  modulation  of  the  surface  wave  magnitude  can 
be  observed  at  a  one-half-wavelength  period.  The  modulation  is  more  clear  in  the  close-up  view 
for  the  perfect  conductor  shown  in  Fig.  3(a).  The  field  structure  gives  the  appearance  of  bound 
modes  that  are  enhanced  when  they  coincide  with  cavities  in  the  surface  profile.  The  surface  wave 
field  structure  for  the  good  conductor  is  similar,  except  it  is  reduced  in  level  somewhat  by  losses 
in  the  metal.  In  Fig.  3(a)  the  field  magnitude  was  computed  on  a  0.2  x  0.2  unit  grid,  giving  a  better 
representation  of  the  surface  profile. 

We  see  that  when  kh  =  0.5  and  kl  =1.0  the  roughness  alone  is  sufficient  to  support  a  surface 
wave,  as  shown  by  the  perfect  conductor  results.  For  the  good  conductor  model,  the  complex 
dielectric  constant  of  the  metal  is  of  secondary  importance  in  determining  the  surface  wave  prop¬ 
erties,  because  in  this  case  the  rough  surface  scattering  losses  are  evidently  more  important  than 
resistive  losses  for  surface  wave  attenuation.  As  a  consequence,  the  surface  wave  properties  of  the 
perfect  and  good  conductor  models  are  similar.  In  contrast,  when  kh  =  0.1  and  kl=  1 .0,  the  rough¬ 
ness  alone  can  support  only  a  weak  surface  wave.  Scattering  losses  are  also  much  reduced,  and  a 
relatively  strong  surface  wave  can  be  produced  on  the  good  conductor,  supported  by  the  complex 
dielectric  properties  of  the  metal. 

Finally,  another  interesting  phenomenon  is  shown  in  Fig.  3(b).  This  figure  shows  the  real  part 
of  the  field  for  the  perfect  conductor,  corresponding  to  Fig.  2(a).  The  color  scale  now  includes 
negative  and  positive  values,  and  the  zero  level  is  in  the  green  region.  The  real  part  of  the  field 
shows  the  phase  fronts,  which  are  advancing  to  the  left.  The  interesting  feature  here  is  the  pres¬ 
ence  of  a  phase  lag  right  near  the  surface  which  occurs  when  the  roughness  is  relatively  large. 
When  kh  =  0.1,  the  phase  fronts  come  smoothly  down  to  the  surface,  showing  no  phase  lag.  This 
is  an  example  of  a  phenomenon  that  can  be  readily  examined  with  field  visualization,  but  may  be 
more  difficult  to  study  with  more  traditional  theoretical  and  numerical  methods. 

In  summary,  field  visualization  can  be  a  very  useful  method  for  gaining  physical  insight  into 
complex  scattering  phenomena  and  should  prove  to  be  a  valuable  supplement  to  theoretical  inves¬ 
tigations  of  electromagnetic  scattering. 


3.  Perturbation  Theory  for  the  Neumann  Boundary  Condition 

The  presence  of  surface  waves  is  related  to  difficulties  that  occur  with  theoretical  treatments 
of  p-polarized  scattering  from  rough  surfaces.  The  general  area  of  interest  here  is  scattering  from 
conducting  surfaces  with  roughness  that  is  small  in  comparison  with  the  radiation  wavelength. 
This  is  the  regime  where  perturbation  theory  is  usually  assumed  to  apply.  However,  it  is  not  com¬ 
pletely  straightforward  to  develop  and  apply  perturbation  theory  to  these  scattering  problems, 
because  of  the  possible  excitation  of  surface  electromagnetic  waves.  The  difficulty  appears  most 
acutely  for  p-polarized  scattering  from  a  rough,  perfectly  conducting  surface  (the  Neumann 
boundary  condition),  which  for  simplicity  we  consider  to  be  a  one-dimensional  (1-D)  surface.  For 
p-polarized  scattering,  the  electric  vector  is  in  the  plane  of  incidence.  In  the  standard  perturbation 
method,  the  scattering  problem  is  solved  using  a  consistent  expansion  in  the  small  parameter  kh, 
where  k  is  the  radiation  wavenumber  and  h  is  the  rms  surface  height  of  the  roughness.  When  this 
approach  is  applied  to  p-polarized  scattering  from  perfect  conductors,  the  lowest-order  result  is 
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incorrect  near  grazing,  and  divergences  appear4,5  beyond  lowest  order.  The  existence  of  surface 
electromagnetic  waves  leads  to  a  breakdown  in  the  standard  perturbation  approach:  Many  of  the 
higher-order  perturbation  theory  terms  are  not  negligibly  small,  even  when  kh  «  1.  A  correct  the¬ 
oretical  treatment  for  this  problem  requires  a  partial  summation  of  these  higher  order  terms;2,5 
that  is,  a  renormalization  procedure  is  required.  (A  large  literature  exists  on  theoretical  treatments 
of  this  problem,  and  no  attempt  will  be  made  here  to  be  complete.)  Because  of  the  need  for  renor¬ 
malization,  the  complexity  of  perturbation  theory  is  significantly  increased  for  this  problem,  and 
numerical  evaluations  are  much  more  difficult. 

For  this  project,  we  have  specialized  to  the  perfect  conductor  case.  The  goal  has  been  to 
develop  a  perturbation  theory  capable  of  giving  accurate  results  for  kh  up  to  0.5.  Perturbation  the¬ 
ory,  if  taken  beyond  lowest  order,  is  valid  for  s-polarized  scattering  from  perfect  conductors 
(Dirichlet  boundary  condition)  up  to  about  this  level  of  roughness.6  However,  renormalization  is 
not  necessary  for  s-polarized  scattering.  The  goal  here  is  to  reach  approximately  the  same  level  of 
accuracy  for  the  Neumann  problem  where  renormalization  is  required.  The  results  to  be  presented 
show  that  this  goal  has  been  essentially  reached. 

Fortunately,  the  fundamental  renormalized  perturbation  theory  for  electromagnetic  scattering 
from  rough  surfaces  has  been  elegantly  worked  out  by  Maradudin  and  co-workers.2"5  However, 
because  of  the  complexity  of  the  results  for  the  scattering  cross  section,  simplifying  approxima¬ 
tions  have  been  used  when  doing  numerical  calculations.  First,  attention  has  focused  on  cases 
with  kh  «  1,  which  allows  a  less  rigorous  evaluation.  Second,  the  pole  approximation  has  often 
been  used  for  the  Green’s  function  for  surface  electromagnetic  waves,  again  simplifying  the  eval¬ 
uation.  What  distinguishes  the  present  approach  is  that  a  full  numerical  implementation  of  the 
renormalized  perturbation  theory  is  made,  without  approximation  beyond  the  basic  perturbation 
approximation  (with  renormalization).  The  most  challenging  part  of  this  program  is  the  solution 
of  a  nonlinear  integral  equation  for  the  “self-energy”  (also  known  as  the  mass  operator).  A  rigor¬ 
ous  evaluation  of  this  quantity  turns  out  to  be  necessary  to  obtain  accurate  results  for  kh  up  to  0.5. 
Previous  work  by  others  has  relied  on  more  approximate  evaluations  of  this  quantity. 

In  Section  3.1  a  summary  of  the  theoretical  results  for  the  scattering  cross  section  and  the 
self-energy  (or  mass  operator)  will  be  given,  together  with  an  outline  of  the  method  used  for 
numerical  evaluation.  A  consistency  condition  known  as  the  Ward  identity  will  also  be  discussed. 
In  Section  3.2,  the  accuracy  of  renormalized  perturbation  theory  for  the  perfect  conductor  case 
will  be  shown  through  comparisons  with  exact  integral  equation  results. 

3.1  Analytical  results  from  perturbation  theory 

We  specialize  to  scattering  from  a  one-dimensional  (1-D)  surface  subject  to  the  Neumann 
boundary  condition  (p-polarized  scattering  from  a  perfect  conductor).  The  strength  of  the  scatter¬ 
ing  will  be  characterized  in  terms  of  the  bistatic  scattering  cross  section6  a.  When  the  incident  and 
scattering  angles  are  measured  relative  to  the  vertical  axis  in  the  usual  fashion,  then  a  is  related  to 
the  differential  reflection  coefficient  dR/dQs ,  also  commonly  used,7  by 

dR 


where  0j  is  the  angle  of  incidence  and  <  >  denotes  an  ensemble  average.  The  plane  of  incidence  is 
the  x-z  plane,  and  the  rough  surface  profile  is  defined  by  z  =fl.x),  where  the  mean  surface  height  is 
zero  and  where  the  z  axis  is  directed  upward  from  the  metal  into  a  vacuum.  We  denote  the  radiation 
wavenumber  by  k,  where  k  =  (0/c.  Horizontal  components  of  wave  vectors  will  be  denoted  by  k 
with  a  subscript.  In  particular,  the  horizontal  components  of  the  incident  and  scattered  wave 
vectors  are  kt  and  ks,  respectively.  Expressions  will  be  given  first  for  standard  perturbation  theory, 
followed  by  those  for  renormalized  perturbation  theory. 

3.1.1  Standard  perturbation  theory 

For  standard  perturbation  theory,  the  lowest-order  result  for  the  scattering  cross  section  is 
given  by 

<*,"  =lu\h.,ks)W(ks-ki).  (11) 

At  lowest  order,  the  field  is  found  to  first  order  in  kh;  thus,  the  scattered  intensity  and  the  scatter¬ 
ing  cross  section  are  to  second  order  in  kh.  This  is  denoted  by  the  1 1  superscript  on  a,  and  the  sub¬ 
script  s  simply  indicates  standard  perturbation  theory.  Also  in  (11), 

u(kvk2)  =k2-kik2  (12) 


and  W(K)  is  a  Gaussian  roughness  spectrum  given  by 


W(K) 


Ih 2  -K2l2/4 


(13) 


where  p(x)  is  the  surface  correlation  function: 


P(*)  =  </(*!)/(*!  +*)>  =  h2 e  X  /l 


(14) 


The  form  given  by  (11)  does  not  vanish  as  0;  or  0S  approach  grazing,  though  it  is  possible  to 
argue  from  quite  general  grounds  that  a  must  have  this  property.  This  indicates  that  (11)  must  be 
inaccurate  close  to  grazing. 

In  order  to  obtain  an  expression  for  the  scattering  cross  section  that  is  accurate  for  greater  val¬ 
ues  of  kh  than  (11),  it  is  logical  to  carry  perturbation  theory  to  higher  order.  When  the  field  is 
found  to  second  order  in  kh,  we  obtain  a  contribution  to  the  cross  section  that  is  fourth  order  in  kh, 
which  we  denote  Gs22.  (Another  term  at  fourth  order,  Gs13,  is  obtained  from  a  product  of  first-  and 
third-order  contributions  to  the  field,  but  this  term  will  not  be  discussed  in  the  context  of  standard 
perturbation  theory.)  The  result  is 

oo 

c22  =  ^jdklW(ks-kl)W(k1-ki)u(ks,kl)u(kl,ki)(j\0]+J20]),  (15) 
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where 


and 


.[0] 


[0], 


=  u(ks,kl)u(k1,ki)  G  (fcj) 


(16) 


Jl20]  =  u(ks,ki  +  ks-kl)u(ki  +  ks-kvki)G[°\kl)Gl0](ki  +  ks-kl),  (17) 


where  *  denotes  complex  conjugation.  In  these  expressions,  the  Green’s  function 


G 


(18) 


where 


%(k} )  — 


i\  k , 


2  2 
k\  <k 

2  2 
k\  >k 


(19) 


As  noted  previously,1’5  the  factor  in  (16)  given  by  IG^^)!2  =  I Jfc2  -  fcj2!'1  causes  the  integral  in 
(15)  to  diverge.  Similarly,  gives  a  divergence  for  backscatter  (ks  =  -kj).  Clearly,  fundamental 
problems  exist  with  standard  perturbation  theory,  and  a  different  method  is  needed. 

3.1.2  Renormalized  perturbation  theory 

Renormalized  (or  self-energy)  perturbation  theory  is  developed  and  discussed  for  scattering 
from  general  media  in  Refs.  2-4  and  specialized  to  the  1-D  Neumann  rough  surface  in  Ref.  5.  A 
renormalized  version  of  a11  is  given  in  Ref.  5,  and  renormalized  expressions  for  a22  and  a13  can 
be  obtained  from  the  general  theory  presented  there.  This  is  accomplished  by  iteration  of  the 
Bethe-Salpeter  equation,  which  is  given  by  Eq.  (6.1)  in  Ref.  5.  Only  the  results  will  be  given  here. 
It  should  be  mentioned  that  the  meaning  of  the  a  superscripts  is  different  for  renormalized  pertur¬ 
bation  theory  than  for  standard  perturbation  theory.  Since  a  partial  summation  of  higher-order 
terms  has  been  accomplished  in  the  renormalization,  contributions  to  all  orders  of  kh  are  present 
in  each  term  of  the  renormalized  perturbation  series.  Nevertheless,  there  is  a  natural  way  of  desig¬ 
nating  a  new  series  of  terms,  which  corresponds  term  by  term  with  the  standard  perturbation 
series;  thus,  the  same  superscript  notation  is  used. 


The  lowest-order  expression  for  the  bistatic  scattering  cross  section  obtained  with  renormal¬ 
ized  perturbation  theory  is  given  by 


11 

a  = 


Ak2 k 2  2 


~[2] 

G  (C) 


G  (ks) 


(20) 


where  kiz  =  -  otp ,(£,■),  ksz  =  ot0(^),  and 
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(21) 


G[2\k  i)  =  - (2) -  ’ 

a Q{kx)-iMK  \kx) 

where  G[2](ki)  is  an  approximation  to  the  surface  electromagnetic  Green’s  function  G{kx),  and 
where  M^2\k{)  is  an  approximation  to  the  self-energy  (or  mass  operator)  M(kx).  (If  the  superscripts 
are  removed  in  (21),  the  proper  relation  connecting  G(k{)  and  M(kx)  is  obtained.)  In  the 
approximation  being  considered  presently,  the  self-energy  is  given  by 


M(2\k  j)  =  jdk2u\kvk2)Gl2\k2)W{kl-k2).  (22) 


Using  (21)  in  (22)  gives  a  nonlinear  integral  equation  for  M^2\kx).  The  renormalization  of  the  the¬ 
ory  leads  to  the  presence  of  M^Qci)  in  (21),  and  it  is  not  difficult  to  see  that  if  M^2\k{)  is  set  to 
zero,  (20)  reduces  to  the  standard  perturbation  theory  result  given  by  (11).  Note  also  that  the  fac¬ 
tors  kj2ks2  in  (20)  lead  to  acceptable  behavior  of  O1 1  near  grazing. 

The  full  self-energy  M(k{)  satisfies  an  integral  equation,1,5  and  approximations  to  M(kx)  are 
obtained  by  iteration,  with  Af2)(£j)  being  the  lowest-order  term.  The  superscript  here  indicates 
that,  aside  from  the  Green’s  function  which  contains  all  orders  of  kh,  the  approximation  to  the 
self-energy  in  (22)  is  of  second  order  in  kh.  Though  (22)  is  the  lowest-order  approximation  to 
M(kx)  obtained  by  iteration,  it  is  still  complicated  to  evaluate  because  a  nonlinear  integral  equa¬ 
tion  must  be  solved.  Because  of  this  complexity,  the  expression  for  A/(2)(kj)  has  apparently 
always  previously  been  further  approximated  by  replacing  G^(&j)  with  G'(>1(£]),  which  reduces 
(22)  to  an  integral.  The  approximation  used  for  M(kx)  determines  the  extent  of  the  processes  that 
are  included  in  the  renormalization  scheme,  and  it  is  not  clear  from  theoretical  considerations 
alone  what  level  of  approximation  is  required.  Some  insight  can  be  gained  by  displaying  the  pro¬ 
cesses  that  enter  using  Feynman  diagrams,2  but  that  will  not  be  discussed  here.  We  found  by 
numerical  evaluation  that  the  further  approximation  to  M^2\kx)  just  mentioned  was  not  adequate 
for  the  relatively  large  values  of  kh  of  interest  here.  Thus,  we  have  used  (22)  as  it  stands,  which 
evidently  has  not  been  done  previously. 

There  is  also  a  consistency  condition  known  as  the  Ward  identity2  that  involves  M(k{)  and 
quantities  that  enter  into  the  scattering  cross  section  series.  This  will  be  discussed  shortly,  but  it 
can  be  mentioned  here  that  (20)-(22)  are  consistent  with  this  condition.  To  obtain  a  renormalized 
perturbation  theory  for  the  Neumann  boundary  condition  that  is  accurate  up  to  kh  =  0.5,  it  is  nec¬ 
essary  to  extend  the  series  beyond  the  o11  term  to  the  a22  and  a13  terms.  We  first  consider  this 
extension  without  improving  the  approximation  to  M(k{),  which  would  be  required  to  remain  con¬ 
sistent  with  the  Ward  identity. 
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(23) 


The  renormalized  results  for  c22  and  a13  are  as  follows 
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(29) 


Comparison  of  (23)-(26)  with  (15)— (17)  shows  that  o22  reduces  to  os22  when  M^{k{)  is  set  to 
zero,  and  a  similar  reduction  occurs  for  a13.  All  expressions  through  (28)  are  valid  for  general 
spectral  density  functions  for  the  surface  roughness,  assuming  the  surface  statistics  are  Gaussian. 
The  last  term  in  (29),  however,  is  specialized  to  a  Gaussian  spectrum.  More  generally,  this  term 
involves  an  integration  over  k2W{k{),  which  can  be  expressed  in  terms  of  the  mean  square  slope 
of  the  surface.  In  (29)  a  Gaussian  spectrum  was  assumed  for  expressing  the  mean  square  slope. 


3.1.3  Evaluation  of  renormalized  perturbation  theory 

The  principal  challenge  in  evaluating  (20),  (23),  and  (27)  is  to  obtain  the  self-energy  M^2\k{). 
Once  that  is  done,  the  remaining  calculation  is  straightforward.  The  natural  approach  to  solving 
the  nonlinear  integral  equation  (22)  is  to  use  iteration,  assuming  M^ikf)  is  small.  Taking  this 
tack,  one  would  initially  set  hf2\k{)  in  G[2\k{)  to  zero  and  iterate  (22).  Unfortunately,  this 
method  does  not  lead  to  convergence  to  a  solution  for  large  Ikj  I  and  does  not  give  the  correct  scat¬ 
tering  cross  section,  except  for  small  kh.  The  problem  is  that  M^2\k{)  is  not  uniformly  small,  but 
instead  grows  quadratically  with  k}  for  large  Ifcjl.  As  kh  increases,  it  becomes  increasingly  impor¬ 
tant  to  compute  M<2>(kl)  accurately  in  this  asymptotic  region.  It  is  not  difficult  to  show  from  (22) 
that  for  very  large  ikjl,  the  leading-order  solution  for  IM^(kj)l  is  k2h.  It  also  follows  from  (21) 
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and  (22),  that  M^2\k{)  must  be  an  even  function  of  k\.  Thus,  using  a  descending  polynomial  form 
in  Ifcjl  beginning  at  the  k{2  term,  we  find  from  (22)  the  following  asymptotic  solution: 


M 


(2) 


(*j) 


=  ik^h  +  -  if  k2h  + 


8  h 


,2 

k  r 


4h 


1  - 


(kiy 


(30) 


Only  the  last  term  involving  kl  depends  on  the  spectral  density  function  being  Gaussian.  All  other 
terms  are  valid  for  arbitrary  W(K).  Equation  (30)  is  valid  for  large  \k}\/k  but  not  for  small  l/ql/k. 
We  found  that  an  adequate  initial  solution  for  M^2\k{)  is  given  by  (30)  when  Im[M0(2)(ki)]  >  0 
and  zero  otherwise.  The  solution  for  frp\k{)  is  then  found  by  repeated  iteration  of  (22). 


For  the  purpose  of  iteration  and  for  later  use  in  (20),  (23),  and  (27),  it  is  convenient  to  repre¬ 
sent  hP-\k{)  as  a  piecewise  continuous  function.  This  was  done  by  selecting  a  finite  range  of  kj 
determined  by  Ikjl  <  k0  over  which  M^kj)  is  to  be  determined.  This  range  was  selected  so  that 
for  larger  \kx\  the  asymptotic  solution  given  by  (30)  is  adequate.  (Typically,  lq,  =  7k.)  The  region 
Ikjl  <  k0  was  then  divided  into  several  intervals,  and  on  each  interval  values  for  the  real  and  imag¬ 
inary  parts  of  A ^2\k{),  computed  from  (22)  with  the  initial  solution  used  in  (21),  were  fitted  by 
40th-order  Chebyshev  polynomials.  This  piecewise  continuous  representation  for  M(1)(k]),  which 
includes  the  asymptotic  solution  for  \k{\  >  k0,  was  then  used  in  (22)  to  obtain  updated  values  for 
M^2\ki),  which  were  again  fitted  with  Chebyshev  polynomials.  The  iteration  process  was  then 
repeated,  and  convergence  to  machine  precision  was  obtained  in  about  ten  iterations.  The  final 
piecewise  continuous  representation  of  M^k^  was  then  used  in  (20),  (23),  and  (27). 

3.1.4  The  Ward  identity 

When  there  is  no  resistive  energy  loss  (as  is  the  case  for  the  perfect  conductor),  the  Ward 
identity2  is  a  relation  that  connects  the  self-energy  M{k{)  with  another  quantity  called  the  aver¬ 
aged  irreducible  vertex  function.2,5  Energy  conservation  will  be  ensured  in  a  given  approxima¬ 
tion  if  the  Ward  identity  is  satisfied.  The  irreducible  vertex  function  satisfies  an  integral  equation 
and  enters  into  the  Bethe-Salpeter  equation,  also  an  integral  equation.  Both  of  these  integral  equa¬ 
tions  are  iterated  in  generating  the  renormalized  perturbation  series  for  the  scattering  cross  sec¬ 
tion.  The  Ward  identity,  then,  is  a  consistency  condition  between  the  approximation  used  for 
M(kx)  and  the  terms  included  in  the  scattering  cross  section  series.  It  is  not  clear  a  priori  how 
important  it  is  to  satisfy  the  Ward  identity  in  order  to  obtain  an  accurate  approximation  to  the  scat¬ 
tering  cross  section. 

The  Ward  identity  is  given  by2 


p  dkn 

Im[ M(k{)]  =  J  2^Im[G(k2)  U(k2,  kj)]  (31) 

— oo 

where  U(k2,  kj)  is  the  averaged  irreducible  vertex  function.  For  01 1  we  need  only  the  lowest-order 
result  for  U(k2,  k{)  given  by 

U{1\k2,kx)  =  2%u\kvk2)  W(kj-k2)  .  (32) 
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With  the  lowest-order  expressions  for  U(k2,  &i)  and  M(k{),  the  Ward  identity  becomes 


Im[M(2)(/t1)]  =  j^u\kvk2)W(kl-k2)Im[G[2](k2)].  (33) 


It  follows  from  (22)  that  the  Ward  identity  is  indeed  satisfied  for  this  case.  It  is  therefore  consis¬ 
tent  with  the  Ward  identity  to  use  (22)  for  the  self-energy  when  computing  a11.  When  the  pertur¬ 
bation  series  for  the  scattering  cross  section  is  extended  to  the  o22  and  o13  terms,  however,  the  use 
of  (22)  is  no  longer  sufficient  to  ensure  satisfaction  of  the  Ward  identity,  because  a  higher-order 
result  is  used  for  U(k2,  k{).  Because  of  the  complexity  of  computing  M(kx)  to  higher  order,  we 
have  at  first  continued  to  use  M^2\kx)  for  a22  and  a13  as  indicated  in  (23)-(28).  Results  showing 
the  accuracy  of  renormalized  perturbation  theory  when  using  this  approximation  will  be  given  in 
the  next  section.  In  addition,  to  examine  the  possible  benefit  of  maintaining  the  consistency 
between  U(k2,  k{)  and  M(k{),  the  self-energy  has  also  been  extended  to  the  next  order,  as 
described  next. 

When  U(k2,  k{)  is  computed  for  the  cross  section  terms  a22  and  a13,  all  terms  to  fourth  order 
in  kh  are  retained,  again  where  only  contributions  apart  from  the  Green’s  functions  are  counted. 
(This  convention  on  order  number  will  be  assumed  henceforth,  and  further  reference  to  excluding 
the  Green’s  functions  will  be  omitted.)  We  therefore  iterate  the  governing  equation  for  M(kx)  and 
keep  all  terms  to  the  fourth  order.  The  result  is  written  in  the  following  form 

M[%!)  =  M 4(2)(it1)  +  M^\k{)  (34) 

where  M4(2)(A:1)  is  the  second  order  part  and  is  given  by  (22),  except  G^2\kx)  is  upgraded  to 
G^fcj),  to  be  defined  shortly.  The  superscript  convention  is  that  (n)  indicates  a  term  of  nth  order 
in  kh,  while  [n]  indicates  that  all  terms  up  to  nth  order  are  included.  The  fourth  order  term  in  (34) 
is  given  by 

M{A\kx)  =  jdk2Gl4\k2)u(kvk2)W(kl-k2)[l(kvk2) +Vl\kvk2)],  (35) 

where 

I(kvk2)  =  \dk3u(k1,k]  +  k3)u(k]  +  kyk3  +  k2)u(k3  +  k2,k2) 

x  G  [4\it1  +  k3)  G  l4\k3  +  k2)  W(k3)  (36) 

and 

VU(kvk2)  =  (kh)2  [2k2  -  (k1+k2)2 +  klk2(k2]+k2)/k2 +  4kxk2/ (kl)2]  .  (37) 
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(38) 


In  these  equations  the  Green’s  function  is  given  by 


~[4L,  ^ 

G  (k})  = 


[4L 


The  Green’s  function  superscript  indicates  the  order  of  the  self-energy  that  is  in  the  denominator, 
and  the  square  bracket  is  used  because  all  contributions  up  to  the  highest  are  included  in  the  self¬ 
energy.  It  can  be  shown  that  the  Ward  identity  (31)  is  satisfied  when  M^4\k{),  G[4\k]),  and  the 
corresponding  lfi4\k2,  k{)  are  used. 

The  computation  of  M[4\k{  )  proceeds  along  the  same  lines  as  described  for  h^2\k{)  and  will 
not  be  discussed  in  detail.  An  asymptotic  form  for  frfi4\ki)  is  used  in  (38)  to  begin  the  iteration  of 
(34)— (38)  in  conjunction  with  (22),  which  is  modified  to  give  MA(1Hkx)  as  mentioned.  After  the 
solution  for  i)  is  obtained,  (20),  (23),  and  (27)  are  evaluated  with  G^(k])  everywhere 
replaced  by  G[4\k{).  Note  that  in  computing  a22  and  a13  such  that  the  Ward  identity  is  satisfied, 
the  computation  of  a11  is  also  modified. 

3.2  Numerical  results  from  perturbation  theory 

In  this  section  comparisons  with  integral  equation  results  are  used  to  examine  the  accuracy  of 
the  various  perturbation  theory  approximations  discussed  in  Section  3.1.  All  examples  assume 
p-polarized  scattering  from  perfect  conductors,  i.e.,  the  Neumann  boundary  condition.  The  sur¬ 
face  roughness  is  described  by  a  Gaussian  roughness  spectrum  with  kl  =  2.0,  and  kh  is  varied.  The 
angle  of  incidence  is  45°  for  all  examples.  The  integral  equation  scattering  cross  sections  were 
obtained  by  averaging  the  results  for  50  surface  realizations,  and  a  tapered  plane  wave  incident 
field  was  used.  Unless  otherwise  noted,  the  surface  realizations  were  80  X  in  length  with  10  parti¬ 
tions  per  wavelength.  We  begin  with  small  kh,  followed  by  steadily  increasing  values. 

Figure  4  shows  the  results  for  kh  =  0.1.  The  fluctuating  curve  was  obtained  with  the  integral 
equation  method,  while  the  smooth  curve  is  the  perturbation  theory  prediction.  Note  that  the 
bistatic  scattering  cross  section  values  are  given  on  a  logarithmic  (dB)  scale.  The  large  peak  at  a 
scattering  angle  of  45°  is  the  coherent  component  in  the  specular  direction;  the  perturbation 
results  are  for  incoherent  scattering  only.  Figure  4(a)  shows  that  the  lowest-order  standard  pertur¬ 
bation  result  [as  11  given  by  Eq.  (11)]  is  accurate  for  this  example,  except  very  close  to  grazing 
(near  ±90°).  For  plane  wave  scattering  from  an  arbitrarily  long  surface,  which  is  assumed  in  per¬ 
turbation  theory,  the  cross  section  expressed  in  linear  units  should  approach  zero  when  the  scatter¬ 
ing  angle  approaches  ±90°.  The  fact  that  it  does  not  is  one  indication  that  a  more  complete  theory 
is  necessary,  as  mentioned  previously.  The  integral  equation  result  drops  down  near  grazing,  but 
does  not  vanish  at  ±90°  because  the  surface  length  is  just  80  X.  When  the  surface  length  is 
increased,  the  cross  section  level  at  grazing  drops  further. 

Figure  4(b)  gives  the  corresponding  result  for  lowest-order  renormalized  perturbation  theory  [Eq. 
(20)].  The  agreement  with  exact  results  is  again  very  good.  Note  that  with  renormalized 
perturbation  theory,  the  prediction  drops  rapidly  very  close  to  ±90°.  The  integral  equation  result 
does  not  have  sufficient  angular  resolution  in  this  region  to  track  the  a1 1  curve.  The  a"2  result  from 
(23)  is  given  in  Fig.  4(c).  For  the  parameters  of  this  example,  the  a22  curve  in  nearly  20  dB  below 
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Figure  4.  Perturbation  theory  compared  with  integral  equation  results  (fluctuating  curve)  for  a 
perfect  conductor  with  kh  =  0.1,  kl  =  1.0,  and  =  45°:  (a)  (b)  cr11 ,  (c)  <T2. 

Renormalized  perturbation  theory  has  been  used,  except  when  the  subscript  s  indicates 
standard  perturbation  theory. 


the  a11  curve  and  can  be  safely  neglected.  (The  a13  curve  is  even  further  below.)  Note,  however, 
the  small  backscattering  enhancement  peak  in  the  backscattering  direction  (-45°).  The 
enhancement  here  is  due  to  the  coherent  addition  of  reciprocal  double-scattering  paths,  where  the 
intermediate  state  is  dominated  by  a  roughness  induced  surface  wave. 

For  Fig.  5,  kh  has  been  increased  to  0.2,  all  other  parameters  remaining  the  same.  Standard 
perturbation  theory,  Fig.  5(a),  is  now  slightly  high  but  still  quite  good.  Recall  that  with  standard 
perturbation  theory  only  the  lowest  order,  os  ,  gives  finite  results,  so  this  series  cannot  be 
extended.  The  renormalized  a11  remains  in  good  agreement  with  exact  results,  as  shown  in 
Fig.  5(b),  though  adding  the  a22  term  [Fig.  5(c)]  does  give  slight  improvement  in  the  forward 
scattering  region  near  the  specular  direction. 

In  Fig.  6,  kh  has  been  increased  to  0.3.  The  surface  length  was  also  increased  to  160  X,  which 
gives  improved  angular  resolution.  The  result  for  osn  [Fig.  6(a)]  has  more  noticeable  error  than  in 
the  previous  example,  and  at  this  level  of  roughness  the  accuracy  of  the  renormalized  a11  [Fig. 
6(b)]  has  also  degraded.  The  result  for  o11  is  about  as  much  low  as  asn  is  high  when  represented 
on  a  logarithmic  scale,  and  this  pattern  continues  at  higher  levels  of  roughness.  It  is  very  reassur¬ 
ing  that  when  a22  is  included  [Fig.  6(c)]  the  agreement  is  excellent.  This  shows  that  the  renormal¬ 
ization  has  properly  converted  the  infinite  cs22  term  into  a  quantity  that  accurately  corrects  the 
error  in  o11  for  this  example.  The  addition  of  a13  does  not  cause  a  noticeable  change  to  the  result 
shown  in  Fig.  6(c). 

Figures  7  and  8  give  the  results  for  kh  =  0.4  and  0.5,  respectively.  The  error  for  a11  is  about 
the  same  as  that  shown  in  Figs.  7(a)  and  8(a)  for  Gsn,  but  of  opposite  sign,  and  is  not  shown. 
Examination  of  the  (b)  and  (c)  panels  of  these  figures  shows  that  the  goal  of  implementing  a 
renormalized  perturbation  theory  valid  up  to  a  kh  of  0.5  has  been  essentially  reached.  For  all  of 
the  examples  discussed  so  far,  the  self-energy  was  computed  according  to  (22).  Therefore,  as  dis¬ 
cussed  in  Section  3.1.4,  the  Ward  identity  is  not  satisfied  when  a22  and/or  a13  are  included. 

Though  the  agreement  shown  in  Figs.  7(c)  and  8(c)  between  renormalized  perturbation  the¬ 
ory  and  exact  results  is  very  good,  there  is  a  little  room  for  improvement,  somewhat  more  notice¬ 
able  in  Fig.  8(c).  The  perturbation  theory  calculations  were  redone  using  the  fourth-order  self¬ 
energy  [Eq.(34)]  to  see  if  satisfaction  of  the  Ward  identity  would  bring  the  results  into  even  better 
agreement.  The  results  are  shown  in  Fig.  9.  When  kh  =  0.4,  the  agreement  is  quite  good  in  both 
Figs.  7(c)  and  9(a),  and  little  can  be  concluded  from  this  comparison.  When  kh  =  0.5,  the  agree¬ 
ment  with  the  Ward  identity  satisfied  [Fig.  9(b)]  is  improved  in  the  forward  region  near  the  specu¬ 
lar  peak  in  comparison  with  the  agreement  when  it  is  not  satisfied  [Fig.  8(c)],  but  is  not  as  good  in 
the  back  region.  For  these  cases,  there  is  no  real  advantage  in  going  through  the  extra  work  of 
upgrading  the  self-energy  in  order  that  the  Ward  identity  is  satisfied  to  fourth  order.  This  can  be 
viewed  as  a  good  result,  lessening  the  complexity  of  extending  these  methods  to  other  problems, 
such  as  scattering  from  good  conductors  or  from  2-D  surfaces.  It  should  be  recognized,  however, 
that  to  achieve  the  success  with  renormalized  perturbation  theory  shown  here,  the  self-energy  was 
required  at  the  full  level  defined  by  (22).  Further  approximation  to  (22)  in  order  to  avoid  solution 
of  the  nonlinear  integral  equation  did  not  yield  satisfactory  results. 
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Figure  9.  Perturbation  theory  result  for  o]]  +  o22+o13  when  the  Ward  identity  is  satisfied  to  fourth 
order  using  the  self-energy  given  by  Eq.  (34):  (a)  kh  =  0.4,  (b)  kh  =  0.5. 
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4.  Summary 

Research  has  been  conducted  in  two  related  areas  for  p-polarized  scattering  from  conducting 
surfaces:  (1)  the  generation  of  surface  waves,  and  (2)  the  development  of  renormalized  perturba¬ 
tion  theory  for  scattering  from  perfectly  conducting  surfaces. 

To  understand  the  excitation  of  surface  waves  better,  a  numerical  method  was  developed  for 
visualizing  the  electromagnetic  field  near  rough  surfaces,  as  illustrated  in  Figs.  1-3.  The  method 
is  based  on  integral  equation  solutions,  and  no  fundamental  approximations  are  required.  Field 
visualization  should  be  very  useful  for  gaining  physical  insight  into  complex  scattering  phenom¬ 
ena  where  theoretical  approaches  may  be  insufficient.  Thus  this  method  has  the  potential  for 
improving  our  understanding  of  scattering  phenomena  in  some  applied  problems. 

The  existence  of  electromagnetic  surface  waves  is  related  to  an  increase  in  the  complexity  of 
theoretical  treatments  of  p-polarized  scattering.  We  have  specialized  to  scattering  from  perfectly 
conducting  surfaces  with  relatively  small  roughness  and  developed  an  implementation  of  renor¬ 
malized  perturbation  theory  that  extends  the  work  of  previous  investigators.  The  accuracy  of  the 
method  was  compared  with  exact  numerical  results,  showing  that  our  implementation  yields  accu¬ 
rate  results  beyond  the  range  of  previous  efforts.  The  work  in  this  project  is  one  step  in  the  more 
general  task  of  developing  robust  scattering  models  for  2-D  surfaces  with  arbitrary  conductivity. 
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